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Abstract 

We present calculations of the magnetic ground states of Cr trimers in different geometries 
on top of a Au(lll) surface. By using a least square fit method based on a fully relativistic 
embedded-cluster Green's function method first we determined the parameters of a classical vector- 
spin model consisting of second and fourth order interactions. The newly developed method requires 
no symmetry constraints, therefore, it is throughout applicable for small nanoparticles of arbitrary 
geometry. The magnetic ground states were then found by solving the Landau-Lifshitz-Gilbert 
equations. In all considered cases the configurational energy of the Cr trimers is dominated by 
large antiferromagnetic nearest neighbor interactions, whilst biquadratic spin-interactions have 
the second largest contributions to the energy. We find that an equilateral Cr trimer exhibits a 
frustrated 120° Neel type of ground state with a small out-of-plane component of the magnetization 
and we show that the Dzyaloshinsky-Moriya interactions determine the chirality of the magnetic 
ground state. In cases of a linear chain and an isosceles trimer collinear antiferromagnetic ground 
states are obtained with a magnetization lying parallel to the surface. 

PACS numbers: 75.30.Hx, 73.22.-f, 75.30.Gw 
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I. INTRODUCTION 



The development of nanoscale devices based on electron spin requires both a fundamental 
understanding of magnetic interactions and practical solutions to a variety of challenges. 
Deposited clusters are of special interest due to their possible application in miniaturized 
data storage technology. The development of scanning tunneling microscopy (STM) and the 
ability to build clusters with well-controlled structures permit the measurement of various 
effects induced by local interactions within magnetic nanoclusters. Recent STM studies 
have investigated the coupling between the magnetic and electronic degrees of freedom of 
nanoparticles and the conducting substrate for adatoms 1 - 1 ^, dimers 4 - 1 ^ and trimers^. Very 
recently Wahl and coworkers^ have been able to estimate the exchange coupling between 
Co atoms on Cu(OOl) surface by probing the Kondo resonance in terms of low temperature 
scanning tunneling spectroscopy. A large number of theoretical efforts has been focused on 
the description of the Kondo effect of single atoms or small clusters.-^ 1 ^ 1 ^ 1 ^ 

First principles studies of supported clusters are often useful for a clear interpretation of 
experimental results and can help a lot in understanding the underlying physical phenomena. 
Determining the, in general, non-collinear magnetic ground states of finite nanoparticles on 
an ab initio level is clearly a demanding task of computational science. One stream of 
such works is based on a fully unconstrained local spin-density approach (LSDA) imple- 
mented via the full-potential linearized augmented plane-wave (FLAPW) method^ or the 
projector augmented-wave (PAW) method 14 . Unconstrained non-collinear calculations are 
also performed within the atomic sphere approximation (ASA) by using a real-space lin- 
earized muffin-tin orbital (LMTO) method 1 ^ 1 ^ 1 ^ and the Korringa-Kohn-Rostoker (KKR) 
method 1 ^. Other works 1 ^^^ rely on ab initio spin-dynamics in terms of a constrained LSDA 
by means of a fully relativistic KKR method and solving simultaneously the Landau-Lifshitz- 
Gilbert equations for the evolution of the orientations of magnetic moments. Although such 
simulations are very accurate in finding the magnetic ground state of complex systems, they 
are very costly and, in practice, require a massively parallel computer architecture. 

Multiscale approaches based on a first principles evaluation of model parameters are very 
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useful to study both the ground state and the dynamics of spin-systems. In Refs. 
the torque method— was employed to calculate isotropic exchange interactions, and then 
Monte-Carlo simulations were performed to study temperature dependent magnetism of nan- 
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oclusters. This approach can, in principle, be extended to include relativistic contributions 
to the exchange interactions^. Nevertheless, because of the low (or even missing) symmetry 
of nanoparticles the determination of the exchange coupling and on-site anisotropy matri- 
ces becomes quite complicated. Moreover, as found, e.g., for Mn and Cr monolayers on 
Cu(lll) higher order spin-interactions are needed for an accurate mapping of the energy 
obtained from first principles calculations^. Recently, a fast ab initio approach that makes 
use of a suitable parametrization of the configurational energy of a complex magnetic sys- 
tem, namely, a spin cluster expansion (SCE), has been proposed^I^, but not yet applied 
intensively. 

In this work we introduce a new scenario to construct parameters of a spin-model con- 
taining interactions, in principle, up to arbitrary order. Our method is based on relativistic 
first principles calculations of the energy, whereby a sufficiently large number of states with 
different non-collinear magnetic configurations (orientational states) are considered to enable 
a least square fit of the parameters of the spin-model. In order to determine the magnetic 
ground state of the system we then solve the Landau-Lifshitz-Gilbert equations derived from 
the corresponding spin Hamiltonian. 

The half-filled valence configuration of Cr yields a large magnetic moment and strong an- 
tiferromagnetic inter-atomic bonding leads in turn to magnetic frustration and complex spin 
phenomena. The simplest system exhibiting such properties is a trimer. The non-collinear 
magnetic structure of supported triangular clusters has been first investigated by means of 
a self-consistent vector Anderson modePi. First principles calculations of an equilateral Cr 
trimer supported on a Au(lll) surface have also confirmed a frustrated non-collinear mag- 
netic structure^ 1 ^ and revealed a collinear antiferromagnetic magnetic ground state for a 
linear chain of three Cr atoms^. 

We apply our new method to Cr trimers deposited on a Au(lll) surface in equilateral, 
linear and isosceles geometries. Though these systems are governed by large antiferromag- 
netic nearest neighbor couplings, we intended to trace the role of the relativistic interactions 
in the formation of the magnetic ground state. Prominently, for an equilateral trimer the 
Dzyaloshinsky-Moriya interactions are shown to fix the chirality of the magnetic ground 
state, whereas in cases of linear and isosceles trimers the inter- and on-site anisotropic terms 
lead to an in-plane orientation of the antiferromagnetic ground state. 
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II. THEORETICAL APPROACHES AND COMPUTATIONAL METHODS 



A. The energy of a classical spin-system 

Neglecting intraatomic non-collinearity, the magnetic state of N atoms is described by the 
array {Mj}^^...^, where M, = di = 1) is the magnetic moment of a particular atom 
labeled by z. In a large class of magnetic systems, referred to as 'good moment' systems, the 
longitudinal fluctuations of the moments can also be neglected, i.e., the magnitudes of the 
moments, Mi, can be considered independent of the orientational state, {<?i}i=i,.,.,jv- The 
most general expression of the energy up to second order of the spin- variables can be written 
as 

E{{a l })=E^ + E^{{a t }), (1) 

with 

E {2) {{Si}) = - * 3 v °i + Bi Ki Bi ' 

i^j i 

where the J i3 - = {J,™ } (a, (3 = x,y,z) are the generalized exchange interaction matrices, 
the Kj = {K"^} are the (second-order) on-site anisotropy constant matrices. Within a non- 
relativistic approach the on-site anisotropy constants vanish, as well as the exchange tensor 
takes a simple diagonal form, Jjj = I with I being the unit matrix, thus, the isotropic 
Heisenberg model is acquired. For a transparent physical interpretation the exchange tensor, 
Jjj, can be decomposed into three terms as^ 

Jy = Jyl + Jg + J* , (3) 

where is the isotropic part of the exchange tensor, 

Jy = ^Tr(jy , (4) 
the traceless symmetric anisotropic exchange tensor, Jg is defined as 

where T denotes transpose of a matrix, while the antisymmetric exchange matrix, J^, as 

J5 = \ (J« - J5) • (6) 



5 



The antisymmetric part of the intersite exchange interaction can then be represented as, 

diXfjOj = Dij{Bi xBj) , (7) 

which is the well-known relativistic Dzyaloshinsky-Moriya (DM) interaction^^, with the 
vector Dij defined as 

D ij = 2 ( J if ~ ■'•!) ' D iJ = 2 ~ ' °*3 = 2 ^ ~ ' ^ 
The asymmetric exchange interactions induced by the spin-orbit coupling have been shown 
to have crucial consequences to the magnetic ground state in thin films.— >2» For transition 
metal clusters such effects are expected to be even more important due to their reduced 
symmetry. 

Unlike most of the thin films with uniaxial or biaxial symmetry, in case of finite clusters 
the structure of the on-site anisotropy matrices can not, in general, be predicted 'a priori', 
i.e., from symmetry principles. The on-site anisotropy can, at best, be characterized by 
diagonalizing the matrix Kj, 

# i K i # i = Y l Kifr-$) 2 , (9) 

A 

where Kf and the unit vectors (A = 1, 2, 3) are the eigenvalues and corresponding eigen- 
vectors of Kj. Clearly, the easy axis is associated by the eigenvector that refers to the 
minimum value of K*. Note that the matrix Kj can be chosen to be symmetric, therefore, 
the eigenvectors are pairwise normal to each other. Obviously, the symmetric anisotropic 
exchange interaction, see Eq. (jSJ), can be decomposed in a similar way, 

Si jg s j = 4 X • e t)(^- • e t) > ( 10 ) 

A 

with J^' x and being the eigenvalues and eigenvectors of the matrix Jg, respectively. 

The second order approximation, Eq. ([2]), is, however, not always sufficient to describe the 
energy of a magnetic system.— As will be shown adding a term, , to Eq. ([1]) corresponding 
to the fourth order spin-interactions considerably improves the quality of the mapping of 
the energy from first principles calculations to the spin-model. In order to keep our model 
tractable, we extended Eq. (jTJ) only by SU(2) invariant fourth order terms, 

e {a \{s 1 })= Qw& ■ °iWk ■ Si) . (n) 

(i<j,k<l) 
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It is easy to see that in case of three atoms the above sum consists of only six different 
terms, therefore, the following simplified notation can be used, 

Ql2 = Ql212, Ql3 = Ql313, Q23 = Q2323, Q23 = Ql213, Q13 = Q2123, Q12 = Q3132 • 

We determined the parameters, J°j , iff , Qij and Q l - k for different Cr trimers on a Au(lll) 
surface by fitting the energy of the orientational states obtained from first principles calcu- 
lations to Eq. (121) supplemented by the terms, Eq. ( TTTT) . 

B. Evaluation of the parameters for Cr trimers 

In order to calculate the electronic structure of the Cr 3 clusters we applied an Embedded 
Cluster Green's function technique as combined with the Korringa-Kohn-Rostoker method 
(KKR-EC)^. Within the KKR-EC the matrix of the so-called scattering path operator 
(SPO), t c , corresponding to a finite cluster C embedded into a host system can be obtained 
from the following Dyson equation, 

rc(E) = T h (E)[l-(t h \E)-t-\E))r h (E)}- 1 , (12) 

where th(E) and Th{E) denote the single-site scattering matrix and the SPO matrix for the 
pristine host confined to the sites in C, respectively, while tc comprises the single-site scat- 
tering matrices of the embedded atoms. Note that Eq. ( fl2i) accounts for all scattering events 
in the system merging the cluster and the host. Once tq is derived, all quantities of interest 
for a cluster, i.e. the charge and magnetization densities, the spin- and orbital moments as 
well as the exchange interaction energy can be calculated. The electronic structure of the 
host gold surface including three layers of empty spheres to represent the vacuum region was 
calculated in terms of the fully relativistic screened Korringa-Kohn-Rostoker method^ 1 ^.. 
The cluster calculations were then carried out such that the Cr atoms substituted empty 
spheres on top of the surface, whereas no attempts were made to include lattice relaxation 
effects. In Fig. [TJ shown is the geometry of the three Cr trimers considered in the present 
work, namely, forming an equilateral triangle, a linear chain and an isosceles triangle. 

The local spin-density approximation as parametrized by Vosko et alM was applied, the 
effective potentials and fields were treated within the atomic sphere approximation. When 
solving the Kohn-Sham-Dirac equation and also for the multipole expansion of the charge 
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FIG. 1: (Color online) Geometry of the Cr trimers (numbered solid circles) deposited on top of a 
Au(lll) surface layer (patterned circles) considered in the present work: a) equilateral triangle, 
b) linear chain, and c) isosceles triangle. The arrows denote the ground state orientation of the 
magnetic moments of the Cr atoms. 

densities we used a cut-off of £ max = 2. When performing self-consistent calculations for 
the linear chain and the isosceles triangle we fixed the direction of the magnetization on all 
the three Cr atom normal to the surface, whilst for the equilateral trimer we used the 120° 
Neel state indicated by arrows in Fig. [T^, as reference (see Sec. III. A). 

For the calculation of the energy of the orientational states we applied the magnetic force 
theorem^^^- by using the self-consistent potentials determined in the reference state. In 
here, only band-energy differences have to be calculated requiring, however, a high precision 
for the necessary Brillouin zone integrals.— To this end, when evaluating Th{E) we used 
over 3300 fc-points in the irreducible (1/6) segment of the Surface Brillouin zone. 

In order to determine the parameters of our spin-model, we generated a large number 
of random magnetic configurations, {c?™}, n — 1, . . . , N, and calculated the corresponding 
band-energies^, Eg = Eb(a™ , a% , a^) . Introducing a (row-) vector containing all combina- 
tions of the components af a (a = x,y,z) occurring in expressions Eqs. and (jTTT) . 

-y I 'n n n n n n n n 

-<^n — \Pl,x a 2,xi a l,x <J 2,yi <J l,x a 2,zi • • • a l,x a l,xi 

• • • > a h a ti a lj a lj"> ' ( 13 ) 

i ,j=x,y,z 

and a vector of the corresponding parameters of the model, 
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the energy of the nth configuration can simply be written as 

E n = PX_ T n . (15) 

An optimal choice of the parameters should minimize the difference (error) between the 
calculated band-energies, E h n , and the energy related to the spin- model, E n . The square of 
error is defined as 

N 

AE* = J2(E n -E b n ) 2 , (16) 

n=l 

and by substituting Eq. (TIB]) a least square condition leads to the solution, 

P = 53^X B fcx^xj) . (17) 

n \ n / 

The number of considered random configurations can be increased until the parameters 
achieve well converged values. 
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FIG. 2: (Color online) The evolution of the relative error, s = \J AE 2 / (E 2 ), of the fitting procedure 
against the number of considered configurations for the isosceles Cr trimer. Solid (red) and dashed 
(blue) lines refer to the parameter set excluding and including SU(2) invariant bi-quadratic terms 
in the spin-model, respectively. 

The quality of the fit is characterized by the relative error, e = V /AE 2 /(E 2 ), where (E 2 ) 
is the average of E\ over all the configurations. This error as a function of the number of 
configurations, N, looks very similar for all the three clusters. As it is shown in Fig. [2] for a 
Cr trimer forming an isosceles triangle the error of the fit is well stabilized around 0.48 % 
above N ~ 5000 when only the second order spin-interactions are taken into account, see 
Eq. (j5J). The error, however, reduced to 0.14 % when also the SU(2) invariant fourth order 
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terms, Eq. (TTTT) . are considered and only about 2000 configurations were sufficient to get a 
stable error. 
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FIG. 3: (Color online) Convergence of the magnitudes of the Dzyaloshinsky-Moriya vectors, D^, 
against the number of configurations for the case of an isosceles Cr trimer. The labels of the DM 
vectors refer to the Cr atoms as numbered in Fig. [lb. 

The configurational energy of the Cr trimers is dominated by quite large antiferromagnetic 
isotropic exchange interactions, ~ 100 — 150 meV. In general, we obtained by two order 
smaller DM interactions, ~ 0.5 — 2.0 meV, whereas the typical range of the anisotropic 
symmetric exchange interactions and the on-site anisotropy constants was about 0.1 meV or 
even less. This trend can be understood in terms of a perturbation treatment with respect to 
the spin-orbit coupling parameter, £, since the DM interactions turn to be proportional with 
£, whereas the anisotropy terms appear (at best) in second order of £. Obviously, a required 
relative accuracy for the small interaction parameters can be achieved at a much larger 
number of configurations than for the total configurational energy. This is demonstrated 
in Fig. [3] showing the evolution of the DM interactions for an isosceles Cr trimer. As can 
be inferred from this figure, about 6-7000 configurations are needed to stabilize the values 
of Dij with a relative accuracy of 1 %. In order to reach the same relative accuracy for 
the coefficients with the smallest magnitude, namely, for the in-plane on-site anisotropy 
constants, we had to generate about 10000 random configurations. 
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C. Determination of the magnetic ground state 



Once the parameters of the spin model are fixed, the ground state configuration of the 
system can easily be determined by solving the Landau-Lifshitz-Gilbert (LLG) equations 
for the transversal components of the magnetizations, 

9Si 7 * . . Tjeff <*7 fjeff 



ir = -t-r-2*< x H i -tjH^ x ft x H i ) > ( 18 ) 

at 1 + a 2 1 + a 2 \ J 

where 7 and a are the gyromagnetic ratio and the Gilbert damping factor, respectively, 
whereas the effective fields, H^, are defined as 

by using Eqs. (j2J) and (fill) . It should be stressed that in the present context Eq. ( fl8i) is 
merely used as a numerical tool to find the energy minimum in the six-dimensional phase- 
space, {<Tj}, describing the non-collinear configurations of the Cr trimers. The stability and 
the speed of the applied numerical procedure, a fourth-order Runge-Kutta method, for in- 
tegrating Eq. (fT51) was, therefore, optimized by adjusting the phenomenological parameters, 
7 and a. 



III. RESULTS AND DISCUSSIONS 



A. Equilateral trimer 

We first investigated a Cr trimer forming an equilateral triangle on top of Au(lll) as 
shown in Fig. [1^,. Since our previous first principles spin-dynamics calculations^ resulted 
in a 120° Neel type of ground state, see Fig. [T^t, in here we used this configuration as a 
reference state to determine effective potentials and fields self-consistently. Reassuringly, 
however, the calculated spin magnetic moments of the Cr atoms, 4.4 fig, were proved to 
be practically independent from the magnetic configuration of the trimer. Since in this 
case both the substrate and the trimer exhibit a c% v point group symmetry, the exchange 
interaction matrices, J 12, J13 and J23, as well as the on-site anisotropy matrices, K 1; K 2 
and K3, are linked in terms of appropriate similarity transformations. The number of 
independent parameters of the model is, therefore, considerably reduced, e.g., the isotropic 
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exchange parameters become identical and the on-site anisotropy matrix corresponding to 
the atom labeled by 2 in Fig. [T^i is of the form, 



K 



XX 











K 2 



K mi K, 



(20) 




zz 



As a consistency check of our fitting process, the obtained parameters satisfied all the rela- 
tionships dictated by the symmetry of the system. 

The dominant parameters determining the ground state of this Cr trimer are the isotropic 
exchange interactions, Jyi — J\3 — J23 — 144.9 meV, and the DM interactions with magni- 
tudes, D12 = D ri = D23 = 1.78 meV. As mentioned earlier, the on-site anisotropy terms are 
much smaller in magnitude, e.g. K xx = —0.09 meV in Eq. (1201 . The DM vectors, visualized 
in Fig. HI point towards a common point lying above the geometrical center of the Cr trimer. 
This is the consequence of Moriya's second rule for the DM vectors^: if a mirror plane of 
the system is bisecting the line between a pair of sites then the respective DM vector lies 
in the mirror plane. Noteworthy, the coefficients of the biquadratic spin-interactions are as 
follows: Q 12 = Q13 = Q23 = -4.42 meV and Ql 3 = Q\ 3 = Q\ 2 = 7.06 meV. 



FIG. 4: (Color online) Schematic view of the DM vectors for an equilateral Cr trimer. 

By solving the LLG equations as described in Section II. C with the above parameters we 
indeed obtained the ground state indicated in Fig. [Da,, namely, a state which is very close to 
an in-plane 120° Neel state with almost negligible out-of-plane components of the magnetic 
moments. Quite obviously, an equivalent ground state can be generated from this state by 
simultaneously reversing the directions of all the magnetic moments. 

Considering only isotropic exchange interactions, a frustration induced by the geometry 
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FIG. 5: (Color online) Two typical ground state configurations of an equilateral Cr trimer in 
absence of DM interactions. The two configurations refer to different chiralities: a) k z = — 1 and 
b) k z = 1, see Eq. (12TJ) . 

of the Cr trimer leads to an eight-fold degenerate non-collinear ground state. These states 
can be divided into two classes as indicated in Fig. [5J One class consists of the configuration 
depicted in Fig. [5^ and the corresponding one with reversed directions (two configurations), 
the second class consists of that in Fig. [5b and those generated from this state via c% v 
symmetry transformations and time reversal (six configurations in total). Defining the 
chirality vector of the system as 

where the summation runs over the three directed bonds, (12, 23 and 31), forming the 
triangle, the two classes can be assigned to chiralities k z = — 1 and k z — 1, respectively. 
Note that for in-plane configurations the vector k is normal to the plane of the triangle. 

Recalling Eq. ([7]), when switching on the DM interactions the degeneracy of the ground 
state will evidently be lifted according to the chirality, k z . In the present case of in-plane 
magnetization the contribution of the DM interaction to the energy can simply be expressed 
as 

E DM = ^-D z k z , (22) 

where D z denotes the z component of any DM vector. Since according to our calculations 
D z = 0.97 meV, the states with k z = —1 depicted in Fig. [5^ become the ground state of the 
system, while the states with k z = 1, see Fig. [5b, are higher in energy by AE = 5.04 meV. 
Thus the ground state we found by solving the LLG equations is a simultaneous consequence 
of antiferromagnetic exchange interactions and the DM interactions as shown in Fig. |H 
Finally in this section, a note has to be added concerning the reference state for the 



fitting procedure described in Sec. II. B. As discussed in quite some details in Ref. 
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when 



choosing a normal-to-plane ferromagnetic reference state an erroneous ground state, Fig. [5b, 
was obtained. The very reason of this result is that in this case the orientations of the fitted 
DM vectors differ from those depicted in Fig. m namely, yielding D z < 0. Clearly from 
Eq. (1221) . the energy of the states related to k z — 1 are lowered with respect to those 
related to k z = —1. This observation clearly indicates that, in particular, for systems with 
metastable states close to the ground state one has to be very careful when choosing the 
reference state serving as basis for subsequent magnetic force theorem calculations. 

B. Linear trimer 

A linear chain of three Cr atoms on top of the Au(lll) surface has been considered in 
the geometry as shown in Fig. [lb. As can be seen from this figure this system has only 
a mirror plane normal to the surface and bisecting the chain. The calculated magnetic 
moments are only slightly different from those for the equilateral triangle: 4.45 fiB at the 
edges of the chain and 4.47 /i# at the central atom. Quite clearly from Tabled! the mirror 
symmetry imposes some relationship between the parameters, e.g., Ju=J23 or J13 has only 
diagonal elements etc. Similar to the equilateral trimer there are large antiferromagnetic 
isotropic exchange interactions between the nearest neighbors, whilst the edge atoms are 
coupled ferromagnetically. The magnitudes of the nearest neighbor exchange interactions 
are very similar to those obtained in terms of a real space LMTO method for Cr dimers^. 
Noticeably, also in this case quite large biquadratic terms of type were needed to obtain 
a sufficiently good fit of the band-energy. 

D 2 1 A ° 23 

\ t / *, 

D 3 1 

FIG. 6: (Color online) Schematic view of the DM vectors for the linear Cr trimer. 

Although interactions of relativistic origin have quite minor contributions to the energy, 
it is instructive to discuss them in some detail. The DM vectors shown schematically in 
Fig. [6] are subjects to symmetry conditions: Di 3 lies in the mirror plane, while D 2 ± and D 2 3 
are mirror images of each other as axial vectors. As in the case of the equilateral triangle, 

14 



Jl2 


7 5 

12 


D21 


Ki 


99.59 


0.157 -0.006 -0.014 
-0.006 0.005 0.063 
-0.014 0.063 -0.162 


0.056 
0.487 
-0.578 


-0.092 -0.007 -0.128 
-0.007 -0.141 0.013 
-0.128 0.013 0.233 




13 


D31 


K 2 


-17.85 


0.022 0.000 0.000 
0.000 0.009 0.000 
0.000 0.000 -0.031 


0.000 
-0.262 
0.068 


0.018 0.000 0.000 
0.000 -0.066 0.034 
0.000 0.034 0.048 


J23 


T 5 

O 90 


^23 


K 3 


99.59 


0.157 0.006 0.014 
0.006 0.005 0.063 
0.014 0.063 -0.162 


0.056 
-0.487 
0.578 


-0.092 0.007 0.128 
0.007 -0.141 0.013 
0.128 0.013 0.233 



Q12 


Q23 


Ql3 


Q23 


Q\3 


Q12 


-0.078 


-0.078 


0.063 


-2.449 


5.810 


-2.449 



TABLE I: Calculated isotropic exchange coupling parameters, Jij, symmetric anisotropic exchange 
tensors, Jfj, DM vectors, Dij, on-site anisotropy matrices, Kj, and SU(2) invariant biquadratic 
coupling parameters, Qij and for a linear Cr trimer. All data are given in units of meV. 

the on-site anisotropy terms and the symmetric anisotropic exchange interactions have the 
smallest contributions to the energy. The on-site anisotropy matrix related to site no. 2 
located in the mirror plane has the structure in Eq. fl20|) . whereas there is no regularity 
in the matrixelements of Ki and K3, except they are related to each other via reflection: 
x' = —x, y' — y, z' — z. 

As expected just by considering isotropic exchange interactions, the solution of the 
LLG equation with sufficiently large damping led to an antiferromagnetic ground state, 
see Fig. [T]d. The direction of the magnetic moments are parallel to the (110) axes which 
is consistent with symmetry considerations, namely, the easy axis of a collinear magnetic 
system with a mirror plane should lie either parallel or normal to the mirror plane^. It is 
important to note that the ambiguity of the reference state mentioned in the previous sec- 
tion does not effect the ground state of the linear chain since the DM interactions evidently 
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vanish in a collinear magnetic state. 
C. Isosceles trimer 

The third type of Cr trimer we considered is an isosceles triangle as depicted in Fig. [lb. 
Apparently, the trimer has a single mirror plane which, however, doesn't coincide with the 
those of the surface layer. Therefore the system has no point-group symmetry in this case. 
Our calculations resulted in magnetic moments of 4.45 \xb for the Cr atoms no. 1 and 3 
and 4.46 /xg for the Cr atom no. 2. As can be inferred from Table HT1 the nearest neighbor 
isotropic exchange interactions are almost symmetric, J 12 — J23, and, as for the linear 
chain, the second nearest neighbor isotropic exchange interaction is weakly ferromagnetic. 
The above data indicate that, similar to many transition metal systems, the formation of 
local moments depends mainly on the nearest neighbor environment of the atoms rather 
than on long-range interactions or the global symmetry of the system. 
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-0.094 0.032 -0.039 
-0.058 -0.039 0.090 


0.610 
-0.335 
1.107 


-0.084 0.024 0.035 
0.024 -0.069 0.006 
0.035 0.006 0.153 



Q12 


Q23 


Ql3 


Q23 


Ql3 


Q12 


1.227 


1.555 


0.271 


-0.640 


3.966 


-0.080 



TABLE II: The same as Table [I] for the isosceles triangle. 
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1 D 31 3 

FIG. 7: (Color online) Schematic view of the DM vectors for the isosceles Cr trimer. 

The second largest contribution to the energy, namely, the biquadratic interactions clearly 
reflect the absence of global symmetry, Qi 2 ^ Q23 and Ql 3 7^ Q\ 2 . This asymmetry is 
more striking in the case of the DM vectors, see also Fig. [7J The good convergence of the 
parameters against the number of configurations seen in Fig. [3] indicates that the large 
asymmetry of the DM vectors indeed stems from the geometry of the system and not from 
an error of the fitting procedure. The obtained antiferromagnetic ground state is very similar 
to that of the linear chain, see Fig. (TJ however, due to the absence of (mirror) symmetry the 
direction of the moments is now slightly out of the line connecting sites no. 1 and 3. 

IV. SUMMARY AND CONCLUSIONS 

We developed a novel method in order to map the energy of supported magnetic nanopar- 
ticles obtained from first principles calculations onto a classical spin Hamiltonian. As a first 
application we determined the spin-interactions for three different Cr trimers deposited on 
a Au(lll) surface. First we calculated the electronic structure of the Cr trimers by means 
of a fully relativistic Green's function embedding method. We obtained magnetic moments 
of the Cr atoms in very good agreement with x-ray magnetic circular dichroism (XMCD) 
measurements^ and with other first principles calculations^ 1 ^. The relativistic treatment 
of the electronic structure was inevitably necessary to properly account for spin-orbit cou- 
pling giving rise to tensorial exchange interactions and magnetic anisotropies influencing the 
formation of non-collinear ground states as shown in case of the equilateral trimer. 

In terms of a least square fit procedure, the most general second-order spin-interactions as 
well as SU(2) invariant fourth-order terms were then fitted serving the best approximation 
to the energies of a large number of random magnetic configurations of the trimers. We 
have shown that the inclusion of fourth-order terms into the spin-model largely enhanced 
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the accuracy of the mapping. A particular advantage of the least square fit applied in this 
work is that it is universally applicable as it does not rely on any symmetry restrictions on 
the model. Moreover, the spin Hamiltonian can be extended to an arbitrary order of the 
spin-interactions. 

The magnetic ground-state of the trimers were found as the solution of the Landau- 
Lifshitz-Gilbert equations. In case of an equilateral Cr trimer we explored that the DM 
interactions lifted the degeneracy of the SU(2) invariant 120° Neel states with different chi- 
rality. On the contrary, for the linear and the isosceles Cr trimers we obtained collinear an- 
tiferromagnetic ground states. An issue of choosing the reference state inherent to methods 
based on the magnetic force theorem was, however, addressed in context to the equilateral 
Cr trimer. This freedom of the method might cause an ambiguity in determining the mag- 
netic ground state of systems exhibiting metastable states close to the ground state. To 
overcome this problem we proposed to use the 'true' ground state obtained from ab-initio 
spin dynamics calculations 2 ^ as reference, since the corresponding spin-model proved to be 
consistent with the 'parent' ground state. 

The present method can be regarded as a very accurate tool in finding the magnetic 
ground state of small supported clusters providing also a clear insight into the role of different 
interactions on the formation of the magnetic ground state. As a prospect for the future, 
the LLG equations will be used to study low-energy spin-excitations of nanoparticles and 
the method can also be extended to include thermal spin-fluctuations.— 1 ^ 
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